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Amplitude  Scaling  of 
Asymmetry-Induced  Transport 

D.L.  Eggleston  and  B.  Carrillo 


Occidental  College,  Physics  Department,  Los  Angeles,  CA  90041 


Abstract.  Our  initial  experiments  on  asymmetry-induced  transport  in  non- neutral 
plasmas  found  the  radial  particle  flux  at  small  radii  to  be  proportional  to  (j>^,  where 
4>a  is  the  applied  asymmetry  amplitude.  Other  researchers,  however,  using  the  global 
expansion  rate  as  a measure  of  the  transport,  have  observed  a (j>\  scaling  when  the 
rigidity  (the  ratio  of  the  axial  bounce  frequency  to  the  azimuthal  rotation  frequency) 
is  in  the  range  one  to  ten.  In  an  effort  to  resolve  this  discrepancy,  we  have  extended 
our  measurements  to  different  radii  and  asymmetry  frequencies.  Although  the  results 
to  date  are  generally  in  agreement  with  those  previously  reported  [(j)\  scaling  at  low 
asymmetry  amplitudes  falling  off  to  a weaker  scaling  at  higher  amplitudes),  we  have 
observed  some  cases  where  the  low  amplitude  scaling  is  closer  to  4i\.  Both  the  and 
(f)\  cases,  however,  have  rigidities  less  than  ten.  Instead,  we  find  that  the  4>\  cases 
are  characterized  by  an  induced  flux  that  is  comparable  in  magnitude  but  opposite 
in  sign  to  the  background  flux.  This  suggests  that  the  mixing  of  applied  and  back- 
ground asymmetries  plays  an  important  role  in  determining  the  amplitude  scaling  of 
this  transport. 


INTRODUCTION 

Malmberg-Penning  traps  are  now  being  used  in  a variety  of  experiments.  Un- 
derstanding the  properties  of  these  traps  is  thus  of  considerable  practical  as  well 
as  fundamental  interest.  It  has  long  been  known  that  confinement  in  Malmberg- 
Penning  traps  is  limited  by  the  presence  of  asymmetric  electric  and  magnetic  fields. 
Early  confinement  studies  found  that  at  low  neutral  pressures  the  confinement  time 
was  much  less  than  expected  from  transport  due  to  electron-neutral  collisions  [1]. 
It  was  suggested  that  this  anomalous  transport  was  due  to  the  presence  of  electric 
or  magnetic  fields  that  break  the  cylindrical  symmetry  of  the  trap.  The  presence  of 
such  asymmetries  would  produce  a radial  component  to  the  E x B drift  that  would 
lead  to  particle  loss  to  the  walls  of  the  trap.  This  suggestion  led  to  a number  of 
experiments  [2-8]  employing  applied  asymmetries  in  order  to  study  the  transport 
in  a controlled  manner.  Most  of  the  experiments  have  used  electric  asymmetries 
since  these  are  easily  applied  and  manipulated  using  the  sectored  wall  portions  of 
the  confinement  region  of  the  trap. 
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A basic  issue  of  asymmetry-induced  transport  is  the  scaling  of  the  transport  with 
asymmetry  amplitude.  Current  theory  [9]  predicts  two  transport  regimes  depending 
on  the  amplitude  of  the  asymmet.ry  in  the  plasma  c/i.  For  smaller  amplitudes  (the 
plateau  regime),  the  radial  flux  is  proportional  to  the  square  of  the  asymmetry 
potential  For  larger  amplitudes  (the  banana  regime),  the  flux  scales  like  the 
square  root  of  the  asymmetry  potential 

Experiments  on  amplitude  scaling  have  used  a variety  of  measures  to  characterize 
the  transport  and  have  applied  the  asymmetries  in  different  ways.  Nott.e  and  Fajans 
[3]  applied  a switched  DC  voltage  to  a 50  degree  wall  patch  in  the  central  part  of 
the  confinement  region.  Following  earlier  work,  they  measured  the  time  t,„  for  the 
central  density  to  decrease  by  one-half  and  found  that  this  confinement  time  scaled 
like  with  1.7  < m < 2.1.  This  scaling  roughly  agrees  with  plateau-regime 
scaling,  although  the  voltages  used  in  this  experiment  (up  to  40  volts)  would  seem 
to  be  much  too  high  to  satisfy  the  plateau  regime  requirement. 

More  recently,  Kriasel  and  Driscoll  [6-8]  have  employed  switched  DC  voltages 
applied  at  the  end  of  the  plasma  column  to  study  this  transport  over  a wide  range  of 
parameters.  For  most  of  these  experiments  they  used  the  rate  of  change  of  the  mean- 
square-radius  Ai/  (with  the  background  rate  subtracted  off)  to  characterize  the 
transport.  They  found  two  transport  regimes  determined  by  the  experimental  value 
of  the  rigidity,  which  is  the  ratio  of  the  axial  bounce  frequency  to  the  azimuthal 
rotation  frequency.  When  the  rigidity  was  between  one  and  ten,  the  expansion  rate 
At*  scaled  like  <f)\,  whereas  higher  values  of  rigidity  gave  a (j>\  scaling. 

Our  earlier  studies  [5]  in  a modified  trap  (described  below)  found  that  the  time 
rate  of  change  of  the  plasma  density  dn/dt  (measured  at  one  radius)  scaled  like 
for  small  asymmetry  amplitude  and  like  for  larger  amplitudes.  The  radial  flux 
r was  found  to  have  a similar  scaling.  The  rigidity  for  these  measurements  was 
small  (approximately  0.2),  and  our  <j)l  scaling  thus  seemed  to  contradict  the  result 
of  Kriesel  and  Driscoll.  The  expansion  rate  Au,  however,  is  related  to  a radial 
integral  [7]  of  the  particle  flux  F,  so  the  contradiction  could  not  be  established  by 
measurements  of  F at  a single  radius.  We  have  thus  expanded  our  range  of  measure- 
ments to  resolve  this  issue.  Although  the  results  are  generally  in  agreement  with 
those  previously  reported,  we  have  observed  some  cases  where  the  low  amplitude 
scaling  is  closer  to  Both  the  <f>l  and  cases,  however,  have  rigidities  less  than 
ten.  Instead,  we  find  that  the  (j)l  cases  are  characterized  by  an  induced  flux  that  is 
comparable  in  magnitiidc  but  ojjposite  in  sign  to  the  background  flux. 


EXPERIMENTAL  RESULTS 

Otir  experiments  are  performed  in  the  modified  Malmberg-Pcnning  trap  shown 
in  Figure  1.  The  usual  plasma  column  is  replaced  by  a biased  wire  running  along 
the  axis  of  the  trap.  Electrons  injected  into  this  device  have  the  same  dynamical 
motions  as  those  in  a normal  non-neutral  plasma  (i.e.  axial  bounce  and  azimuthal 
drift  motions),  but  collective  contributions  to  the  asymmetry  potential  are  largely 
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FIGURE  1.  Schematic  of  the  Occidental  Trap.  The  plasma  is  replaced  by  a biased  wire  that 
maintains  the  basic  dynamical  motions  of  the  injected  electrons.  Forty  wall  sectors  allow  for  the 
application  of  asymmetries  consisting  of  essentially  one  Fourier  mode. 


eliminated  since  the  lower  density  (10^  cm“^)  and  higher  temperature  (4  eV)  of  the 
electrons  give  a Debye  length  larger  than  the  trap  radius.  The  asymmetry  poten- 
tial within  the  confinement  region  is  thus  essentially  the  vacuum  potential  and  we 
need  not  be  concerned  with  plasma  modifications  of  the  asymmetry  potential.  De- 
spite these  changes  in  the  plasma  parameters,  the  confinement  time  scaling  with  no 
applied  asymmetries  [10]  shows  the  same  {L/B)^  dependence  found  in  higher  den- 
sity experiments  [11,12],  thus  supporting  the  notion  that  the  asymmetry-induced 
transport  is  a single  particle  effect. 

For  the  current  experiments,  forty  wall  sectors  are  employed  to  produce  the 
applied  asymmetry.  By  judiciously  selecting  the  voltage  applied  to  each  sector 
we  can  produce  an  asymmetry  consisting  of  essentially  a single  Fourier  mode,  thus 
eliminating  the  sum  over  modes  in  the  transport  theory  [9]  and  making  for  a simpler 
comparison  between  theory  and  experiment.  We  use  AC  voltages  at  a variable 
frequency  / since  this  gives  us  an  additional  experimental  parameter  that  can 
be  varied  independently  of  other  quantities.  Electrons  injected  into  the  trap  are 
quickly  dispersed  into  an  annular  distribution  [13].  At  the  end  of  an  experimental 
cycle  the  electrons  are  dumped  onto  a phosphor  screen  and  the  resulting  image  is 
digitized  with  a cooled  CCD  camera.  A radial  cut  through  this  image  gives  the 
axially-integrated  density  profile  n(r,  t)  of  the  electrons.  For  this  study  it  is  then 
useful  to  numerically  integrate  the  density  profile  to  obtain  the  quantity  g{r,t): 

g{r,  t)  = — / rdr  ■ n(r,  t) 
r Jo 

The  time  derivative  of  g{r,t)  then  gives  the  radial  particle  flux  F(r,  t)  = dg/dt. 

Figure  2 shows  a typical  set  of  data.  The  quantity  g{r,t)  is  plotted  vs  time  for 
various  asymmetry  amplitudes.  A line  is  then  fit  to  the  initial  slope  of  these  curves 
to  obtain  the  radial  flux.  Note  that,  at  the  higher  amplitudes  and  later  times,  the 
slope  changes  as  the  initial  plasma  is  modified  by  the  transport.  Plotting  the  data 
in  this  way  allows  us  to  note  and  avoid  these  saturated  cases.  Also  note  that  the 
flux  is  not  zero  for  zero  applied  asymmetry,  i.e.  there  is  a nonzero  background 
flux  Fq.  Finally,  note  that  because  we  have  an  annular  plasma  both  positive  and 
negative  fluxes  are  observed. 
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FIGURE  2.  A typical  data  set  sliowing  the  computed  quantity  g{r,t)  versus  time  with  iisyni- 
metry  amplitude  as  a parameter.  The  slope  of  these  curves  gives  the  radial  flux  F = dgfdt. 


Figure  3 shows  the  magnitude  of  the  net  flux  AF  = IFinriuccd  “ Fo|  vs  asymmetry 
amplitude  with  radial  position  as  a parameter.  Tlie  experimental  conditions  arc 
similar  to  those  of  reference  [5],  figure  3 (center  wire  bias  = -80V,  B = 365G,  / 
= IMHz).  The  log-log  plot  allows  us  to  determine  the  scaling  exponent  m,  whore 
AF  oc  0”.  Typical  error  bars  are  shown,  and  lines  of  slope  one  and  two  arc  drawn 
for  comparison  to  the  data.  The  results  are  similar  to  those  reported  in  reference 
[5]:  the  net  flux  scales  like  0^  for  low  amplitude  asymmetries,  in  agreement  with 
plateau-regime  theory.  At  higher  amplitudes,  the  scaling  exponent  falls  off  to  a 
smaller  value  in  the  range  0.5  to  1.4.  This  may  indicate  a transition  to  bananei- 
regime  scaling,  but  the  data  is  not  sufficient  to  support  this  conclusion.  Note  that 
there  is  no  dramat  ic  dependence  of  m on  radius,  so  if  we  were  to  calculate  the 
expansion  rate  Au  as  in  reference  [8]  the  low  amplitude  scaling  would  still  be  0„. 

Further  exploration  of  our  parameter  space  has  found  some  cases  showing  scaling 
closer  to  <p\  and  some  of  these  cases  are  shown  in  Figure  4.  For  these  parameters 
(as  above  except  / — 0.35MHz),  the  0‘  cases  occur  at  smaller  radii  while  (j>l  cases 
occur  at  larger  radii,  so,  again,  a calculation  of  expansion  rate  would  not  yield  a 
<j>l  scaling. 

Under  what  conditions  do  we  obtain  fluxes  with  a cpl  scaling?  In  Figure  5 we 
plot  the  scaling  exponent  m versus  the  rigidity  for  the  cases  wo  have  considered. 
Note  that  all  these  cases  have  ridigity  less  than  ten  (the  borderline  value  for  the 
Kriesel/Driscoll  regimes),  and,  more  importantly,  that  the  same  value  of  rigidity 
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FIGURE  3.  Log-log  plot  of  the  net  asymmetry-induced  flux  AF  vs.  the  applied  asymmetry 
amplitude  at  the  wall  <j>a.  The  asymmetry  frequency  f is  IMHz.  Lines  of  slope  one  and  two  are 
shown  for  comparison.  For  the  cases  shown  the  low  amplitude  scaling  is  close  to  0^. 
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FIGURE  4.  Log-log  plot  of  AF  vs.  for  f = 0.35MHz.  For  the  lower  radii  the  scaling 
exponent  is  close  to  one. 
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FIGURE  5.  Scaling  exponent  m vs.  rigidity.  No  correlation  is  oliserved  with  this  parameter. 


gives  both  m ^ 1 and  m w 2 cases.  Thus,  in  our  experiment  the  scaling  exponent 
is  not  correlated  with  rigidity. 

We  have  tried  to  find  correlation  with  other  parameters  as  well.  Since  our  exper- 
iments use  a nonzero  asymmetry  frequency  /,  it  might  make  sense  that  the  rigidity 
should  be  modified  from  fb/fn  to  fb/{f  - /r)  since  presumably  this  ratio  measures 
rigidity  relative  to  the  asymmetry.  A plot  of  the  data,  however,  shows  no  corre- 
lation of  the  scaling  exponent  with  this  modified  rigidity  either.  Other  quantities 
that  fail  to  correlate  with  m include  / — /r,  (/  — fa)/ To,  and  Tinduced- 

A re-examination  of  the  original  data  for  the  m = 1 cases  shows  typical  gvst  plots 
like  the  one  shown  in  Figure  6.  For  these  cases,  the  background  flux  is  comparable 
in  magnitude  to  the  induced  flux.  Following  this  clue,  we  have  plotted  in  Figure  7 
the  scaling  exponent  vs  the  ratio  of  the  background  flux  to  the  induced  flux  at  a 
typical  asymmetry  value  of  0.3V.  The  correlation  here  is  good,  with  m = 2 cases 
corresponding  to  smaller  values  of  Fo/Fo  av  and  the  m = 1 cases  corresponding  to 
larger  negative  values. 


DISCUSSION 

The  fact  that  our  m = 1 cases  occur  when  the  background  flux  is  comparable  in 
magnitude  to  the  induced  flux  suggests  a simple  mixing  of  applied  and  background 
asymmetries.  Suppose  that  the  true  flux  scaling  is  F oc  (j)^  and  that  when  both 
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a background  asymmetry  (po  and  an  applied  asymmetry  </>„  are  present  they  add 
directly.  The  total  asymmetry  squared  is  then  (p^  + 2(pa(po  + (pQ  and  U <pa  « (pQ, 
we  obtain  F - To  oc  epa-  Such  a model,  however,  is  not  consistent  with  the  details 
of  the  experimental  data  or  with  the  theory.  At  least  one  experimental  case  where 
the  flux  scales  like  4>\  occurs  in  the  middle  of  our  plasma  with  (p1  cases  at  smaller 
and  larger  radii.  The  above  model  would  then  require  a background  asymmetry 
that  is  peaked  in  the  middle  of  the  plasma  and  it  is  hard  to  imagine  how  such  an 
asymmetry  could  be  produced.  This  direct  addition  model  is  also  inconsistent  with 
the  theory,  which  says  the  flux  is  given  by  a sum  of  terms,  one  for  each  Fourier 
mode  produced  by  the  asymmetries.  In  our  case,  the  background  asymmetry  is 
static  (i.e.  zero  frequency)  while  the  applied  asymmetry  has  a nonzero  frequency. 
These  will  necessarily  give  different  Fourier  modes  that  should  not  add  directly.  It 
would  seem  that  the  explanation  for  our  observed  correlation  must  be  at  a deeper 
level  than  simple  addition  of  asymmetry  potentials. 

It  is  interesting  to  ask  whether  the  Kriesel/Driscoll  results  can  be  explained 
in  terms  of  our  findings.  From  the  data  in  reference  [7],  it  appears  that  the  <pl 
scaling  does  occur  for  cases  where  the  induced  transport  is  comparable  to  the 
background  transport.  The  <pl  cases  arc  harder  to  judge  because  the  background 
was  so  small  it  was  declared  to  be  zero.  Although  this  would  fit  into  our  scheme 
(i.e.  ro/Vinduced  would  be  zero),  a more  careful  measurement  would  be  required  to 
settle  this  question. 

It  has  been  suggested  [14]  that  the  (pl  scaling  is  due  to  the  large  size  of  the 
asymmetries  employed  in  some  experiments.  In  these  cases  the  small-perturbation 
assumption  of  the  quasi-linear  theory  is  violated.  If  the  applied  asymmetry  is  large 
enough  {ecp  » kT),  particles  will  be  excluded  from  (or  pulled  into)  the  vicinity 
of  the  biased  wall  sector.  The  perturbation  of  the  distribution  function  <5/  would 
thus  be  large  and  saturated  (i.e.  no  longer  dependent  on  the  perturbed  potential 
in  the  plasma  6<p).  Since  the  transport  goes  like  the  product  6f6(p  the  dependence 
would  be  linear  in  the  asymmetry  amplitude.  While  this  model  might  explain 
some  experimental  results  (for  example,  the  fall-off  of  our  scaling  exponent  with 
asymmetry  amplitude),  it  does  not  seem  consistent  with  the  bulk  of  the  data,  either 
in  our  experiments  or  those  of  others.  For  example,  when  we  observe  (pl  scaling,  it 
is  observed  even  at  very  low  asymmetry  amplitudes.  In  contrast,  Notte  and  Fajans 
observed  <pl  scaling  even  at  extremely  high  asymmetry  potentials.  It  is  also  not 
understood  why  this  model  should  depend  on  rigidity,  so  it  cannot  explain  the  two 
scaling  regimes  observed  by  Kriesel  and  Driscoll. 


CONCLUSION 

We  find  that  in  most  cases  the  asymmetry-induced  transport  in  our  experiment 
produces  a net  radial  flux  that  scales  like  the  square  of  the  applied  asymmetry 
amplitude  for  small  asymmetries  and  falls  off  to  a weaker  dependence  (m  = 0.5  - 
1.4)  for  larger  amplitudes.  We  have  found  some  cases  that  give  a low-amplitude 
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scaling  exponent  closer  to  one,  and  these  cases  occur  when  the  induced  flux  is 
comparable  in  magnitude  but  opposite  in  sign  to  the  background  flux.  We  do  not 
yet  understand  the  reason  for  this  correlation. 
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